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The nonlinear and chaotic behaviour of classical field theories is currently 
subject of intensive research [1-3] and, in this respect, it is of great interest to 
investigate the existence and properties of limit cycles, which are inherently 
nonlinear phenomena [4,5]. 

In a previous paper [6] we studied the stability of a scalar infiaton field 
with a symmetric double-well self energy. We showed that the value of the 
inflaton field in the vacuum is a bifurcation parameter which changes the 
phase space structure and that for some functional solutions of the Hubble 
"constant" the system goes to a hmit cycle, i.e. to a periodic orbit. 

In this paper we analyze the properties of this limit cycle by using ana- 
lytical and numerical arguments. We show that the limit cycle is unique and 
stable and give an analytical formula for its period. 

To solve the three major cosmological problems, i.e. the flatness problem, 
the homogeneity problem, and the formation of structure problem, it is gen- 
erally postulated that the universe, at a very early stage after the big bang, 
exhibited a short period of exponential expansion, the so-called inflationary 
phase [7-10]. All the inflationary models assume the existence of a scalar 
held 0, the so-called inflation fleld, with Lagrangian [9,11] 

L = ^d^<i>d^<t> - y (0) (1) 

where the potential V^(0) depends on the type of inflation model considered. 
Here we choose a real fleld but also complex scalar can be used [9]. The 
scalar fleld, if minimally coupled to gravity, satisfles the equation 

□0 = + 3(^)0 -lvV = -^, (2) 
a 0(p 
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where □ is the covariant d'Alembertian operator and a is the cosmological 
scale factor. The parameter G is the gravitational constant {G — with 
Ti — c — 1 and Mp — 1.2 • 10^^ GeV the Plank mass) and H — d/a is the 
Hubble "constant", which in general is a function of time (Hubble function). 
Wc suppose that in the universe there is only the infiaton field, so the Hubble 
function H is related to the energy density of the field by 



Immediately after the onset of inflation, the cosmological scale factor a grows 
exponentially [9]. Thus the term V^0/a^ is generally believed to be negligible 
and, if the inflaton field is sufficiently uniform (i.e. 0^, (V0)^ << V{(j))), we 
end up with a classical nonlinear scalar field theory in one dimension 

+ 3i/(0)0+^ = O, (4) 

where the Hubble function H satisfies the equation 

= ^VW. (5) 

The potential V{(f)) depends on the type of inflation model considered 
and we choose a symmetric double-well potential 

y(0) = ^(0^-.^)^ (6) 

where ±v are the values of the inflaton fleld in the vacuum, i.e. the points 
of minimal energy of the system. 

One of the main difficulties in constructing models with potentials suit- 
able for inflation is that these potentials must be flat enough to allow a 



sufficiently long period of inflation [8,9]. In this respect our model is very 
schematic but it can be seen as a toy model for classical nonlinear dynamics 
with the attractive feature that it emerges form the inflationary cosmology. 
Obviously the complete study of the dynamics of the inflaton may be ad- 
dressed only by a complete quantum field theory approach able to predict 
not only the behaviour of the classical value of the inflaton field but also the 
associated quantum fluctuations. However the quantization of the inflation 
scenario is still an open problem [10] (an interesting stochastic approach can 
be found in [12]). 

In the paper [6] we showed that the inflaton field value in the vacuum v is 
a bifurcation parameter. If f = in the phase space there is only one stable 
fixed point (0 = 0, ^ = 0), which is an attractor. Instead lor v ^ Q there are 
three fixed points: (0 = 0, = 0), which is unstable, and (0 = ±i>, = 0), 
which are stable. 

The Hubble function is determined by solving the equation (5). There 
are four possible continuous solutions : 

H{ct^) = ±^\<t^'-v\ (7) 

but also: 

H{ct>)^±^{ct>'-v% (8) 

where 7 = ^J2'kGX/2> is the friction parameter. The choice of the solution is 
crucial for the dynamical evolution of the system. 

By using the Bendixon criterion [13] (discussed in detail in [6]) we obtain 
that if if(0) = 7|0^ — w^l then the Hubble function does not change sign 
and we do not find periodic orbits. The 4th-order Runge-Kutta numerical 
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integration [14] of the equations of motion shows that for w 7^ the inflation 
field approaches one of its two stable fixed point attractors, and that the 
Hubble function goes to zero with an oscillatory behaviour (see Figure 1). 
Instead, if we choose H{(f)) — 7(0^ — v^) then the Hubble function can change 
sign and the Bcndixon criterion admits for f 7^ the existence of a limit cycle. 
In fact, the numerical calculations plotted in Figure 2 show that a limit cycle 
exists and that the Hubble function oscillates forever. 

Now we want analyze in detail the properties of this limit cycle. The 
equation of motion of the inflaton field with H{(f)) — ^{(fP — v'^) reads 

+ 37(0^ - v'^)4) + A0(02 - = . (9) 

This equation can be written as 

^[0 + 37 _ ^2^^^| ^ ^^^^2 _ ^ Q ^ (^Q) 

at Jo 

and if we put 

FU) =3 ['^(u'^ -v^)du = M'^ -3v^) , GU) = M'^ - v^) , (11) 
^0 

and also uj = (f) + ^F{(f)), we obtain the system 

(f) — uj — 7^(0) 

< 

uj = -XG{(t)) . 

For systems of this kind the Lienard theorem [15,16] states that there is an 
unique and stable limit cycle if the following conditions are satisfied: F{(f)) 
is an odd function and F(0) = only at = and (j) = ±0;; F{(j)) < for 
< < a, F{(j)) > and is increasing for (p > a; G{(j)) is an odd function 
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and (/)G{(f)) > for all (p > a. It is easy to check that the functions 
and defined by (11) satisfy all the conditions of the Lienard theorem 

with a — V. The cubic force tends to reduce any displacement for 

large |0|, whereas the damping F{(/)) is negative at small |0| and positive at 
large |0|. Since small oscillations arc pumped up and large oscillations are 
damped down, it is not surprising that the system tends to Seattle into a 
self-sustained oscillation of some intermediate amplitude. 

Figure 2 and Figure 3 show that both internal and external initial condi- 
tions generate trajectories which approach the limit cycle, so we have also a 
numerical evidence of the stability of the limit cycle. 

Let us consider a typical trajectory of the Lienard system (12). After the 
scaling = Ao; we obtain 

^' = A[^-iF(0)], ^^^^ 

The cubic nuUchne ijj = (7/A)F(0) is the key to understand the motion [5]. 
Suppose that A >> 1 and the initial condition is far from the cubic nuUcline, 
then (13) implies \4>\ ~ 0{X) » 1; hence the velocity is enormous in the 
horizontal direction and tiny in the vertical direction, so trajectories move 
practically horizontally. If the initial condition is above the nuUcline then 
4> > 0, thus the trajectory moves sideways toward the nuUcline. However, 
once the trajectory gets so close that i/j ~ (A/7)F(0) then the trajectory 
crosses the nuUcline vertically and moves slowing along the backside of the 
branch until it reaches the knee and can jump sideways again. The period 
T of the limit cycle is essentially the time required to travel along the two 
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slow branches, since the time spent in the jumps is neghgible for large A. By 
symmetry, the time spent on each branch is the same so we have 

rtB 



T ~ 2 

ItA 



/ dt (14) 

JtA 



where A and B are the initial and final points on the positive slow branch. 
To derive an expression for dt we note that on the slow branches with a good 
approximation ip ~ {^/X)F{(j)) and thus 

Since from (13) dip/dt — —0(0^ — v^), we obtain 

dt ^ -Sl-^l , (16) 

A 

on the slow branches. The slow positive branch begins at (f)A — 2jv/X and 
ends at (pB — ^v/X, hence 

r~2 r^'dt- -6^ T''^ ~6^1n2. (17) 

JtA X J(I>A (p X 



Because 7 = W27rGA/3 we have 



Tc.21n2^^. (18) 

Note that the period is v independent. 

In summary, we have proved the existence and stability of a limit cycle 
in the phase space of a scalar inflaton field (f) with a symmetric double-well 
potential V{(f)) and a friction term in the equation of motion proportional to 
y(0). Then we have obtained an analytical estimation of the period of the 
limit cycle. 
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Figure Captions 



Figure 1: The Hubble function vs time (up) and the phase space trajectory 
of the inflaton field (down); for H{(f)) = 7|0^ — v'^\ with 7 = 1/2, A = 3 and 

V — 1. Initial conditions: = and (j) — 1/2. 

Figure 2: The Hubble function vs time (up) and the phase space trajectory 
of the inflaton field (down); for H{4>) — 7(0^ — i;^) with 7 = 1/2, A = 3 and 

V = 1. Initial conditions: = and = 4. 

Figure 3: The Hubble function vs time (up) and the phase space trajectory 
of the inflaton field (down); for H{(f)) — 7(0^ — v'^) with 7 = 1/2, A = 3 and 

V — 1. Initial conditions: = — 1/2 and = 0. 
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